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ABSTRACT 

Recent advances in supercomputers and highly scalable 
quantum chemistry software render computational 
chemistry methods a viable means of providing chemistry 
data for aerothermal analysis at a specific level of 
confidence. Four examples of first principles quantum 
chemistry calculations will be presented. Study of the 
highly nonequilibrium rotational distribution of a nitrogen 
molecule from the exchange reaction N + N 2 illustrates 
how chemical reactions can influence rotational 
distribution. The reaction C 2 H + H 2 is one example of a 
radical reaction that occurs during hypersonic entry into 
an atmosphere containing methane. A study of the etching 
of a Si surface illustrates our approach to surface 
reactions. A recently developed web accessible database 
and software tool (DDD) that provides the radiation 
profile of diatomic molecules is also described. 

1. INTRODUCTION 

In an aerocapture mission, there is a tradeoff between the 
mass of the propellant needed for orbital insertion and the 
mass of the Thermal Protection System (TPS) needed to 
shield against aerodynamical heating. Thus the design of 
such a mission must include careful aerothermal analysis 
carried out at a specific level of confidence. This in turn 
calls for a corresponding level of confidence in the 
chemistry data used in the analysis. 

The chemistry data needed in aerothermal analysis 
include chemical reaction rates, radiation emission and 
absorption probabilities, and probabilities of surface 
catalytic reactions. There are two sources for these data 
— laboratory measurements, and first principles quantum 
chemistry calculations. Traditionally, experimental 
measurements were the major source of data. However, 
measurements of high temperature radical reactions and 
photoemission probabilities for high-lying rotational and 
vibrational states of a molecule pose experimental 
problems due to the difficulties in the unambiguous 


measurement of a specific reaction rate or molecular 
transition probability at a high temperature. Under such 
circumstances, first principles calculations frequently are 
the only source of reliable data. While first principles 
calculations have been used in the past to provide selected 
data, the approach has been plagued by the lengthy turn- 
around time. Recently, the availability of powerful new 
supercomputers such as Project Columbia at Ames 
Research Center and the development of highly scalable 
chemistry software make quantum chemistry calculations 
more efficient than those of even just a few years ago. 
Thus first principles calculations can now be used as a 
reliable source of chemistry data for aerothermal analysis. 

To illustrate recent advances in modern quantum 
chemistry methods and how they can be used to serve the 
need of mission design for planetary probes, this paper 
presents examples of four recent chemistry calculations. 
The rotational temperature of N 2 is a determining factor of 
its radiation profile and has to be included in the 
calculation of the radiative heat load to a vehicle or probe 
during high-velocity entry into a nitrogen-containing 
atmosphere. Our calculation of the N + N 2 exchange 
reaction is carried out to reconcile the rotational 
temperature of N 2 observed in a recent shock tube 
experiment [1]. The ethynyl radical C 2 H is an important 
intermediate product in the entry shock of a methane 
containing atmosphere and the H 2 + C 2 H reaction 
presented in this paper illustrates current quantum 
chemical capability in calculating high-temperature 
radical reactions. The etching of a silicon surface by 
halogens serves as an example of a computational study 
of surface catalytic reactions. A newly developed web- 
based tool, Dynamic Database for Diatomics (DDD) 
provides the user with the unprecedented capability of 
determining the complete radiation profile of a diatomic 
molecule. 

Even with recent software and hardware advances, the 
amount of chemistry data required in aerothermal analysis 
is still far too numerous for a reasonable turn-around time. 
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In this respect, the recent sensitivity study of 
thermochemical modeling by Bose et al. [2] for Titan 
atmospheric entry points to sensitivity analysis as an 
important tool in selecting a subset of chemical data 
particularly important in reducing the uncertainties in the 
modeling. Sensitivity studies serve as an important guide 
to selecting the needed data before quantum chemistry 
calculations are carried out. 

2. N + N 2 EXCHANGE REACTION 

In a recent shock tube experiment using nitrogen gas, 
Fujita et al. [1] measured the rotational and vibrational 
temperatures using the N 2 second positive and N 2 + first 
negative bands and found significant nonequilibrium 
between the measured rotational and translational 
temperatures. So far, proposed explanations [3,4] of this 
nonequilibrium rotational distribution using N 2 -N 2 
collisions require physically unrealistic features in the 
N 2 -N 2 interaction potential and are likely to be unreliable 
[5]. Thus it is likely that the large A / transitions indicated 
in the shock tube data of Reference [1] do not come from 
N 2 -N 2 collisions. Instead, our study [5] shows that the N 
+ N 2 exchange reaction can lead to large A j transitions at 
the high-temperature regime of the experiment. 

In the atom exchange reaction 

N + N 2 (vj) ^N 2 (v '/) + N, ( 1 ) 

the N-N bond in the original N 2 molecule is broken and a 
new N-N bond is formed. Due to the nature of the N + N 2 
interaction potential [6], the new N 2 molecule may result 
in a highly excited rotational state in comparison with the 
rotational state of the original molecule. 

A theoretical determination of exchange scattering must 
be based on a quantum mechanical description. For 
example, the N-N + resonance charge exchange [7] would 
be underestimated by almost an order of magnitude using 
classical mechanics. The first ab initio potential energy 
surface for the N-N 2 system (denoted as WSHDSP PES) 
that includes variation in the interatomic separation 
distances of the N 2 molecules and the first quantum 
dynamics study of the cumulative reaction probabilities 
and chemical reaction rates of the N + N 2 exchange 
reaction have recently been reported by our group [6]. 

The WSHDSP PES [6] for N-N 2 is constructed from the 
results of high-level quantum electronic structure 
calculations. These calculations are primarily based on the 
same methods [8] that were used earlier to determine the 
ab initio data for the construction of the very accurate 
N-N 2 rigid-rotor PES of Stallcop et al. [9]; long-range 
dispersion force data and nearly 4,000 ab initio data 
points are applied to construct the N-N 2 WSHDSP PES. 
This potential has a number of special features. In 
particular, strong enhancements in the collision cross 


section due to the formation of a short-lived metastable 
N 3 state have been reported [10,1 1]. 

2.1 Rotational Excitation 

State-to-state cross sections for scattering energies at 2.5 
eV are presented in Fig. 1 to illustrate the distribution of 
rotational state excitation. Note that we plot two sets of 
cross sections — one with nuclear spin of the N 3 system 
properly accounted for, and the other without spin. N 3 
obeys the Bose-Einstein statistics. Since the incorporation 
of nuclear spin introduces rapid oscillations in the cross 
sections, the underlying features in the cross sections are 
more easily recognizable using the data without nuclear 
spin. Note, however, the rapidly varying cross sections of 
even versus odd j states in N 2 will be reflected in its 
radiation profile. 



FIG. 1 . State-to-state cross sections for exchange scattering with 
rotational excitation from initial ground rotational and 
vibrational state of N 2 , N + N 2 (v=0, 7=0) — > N 2 (v'=0, j') + N, 
vs. j' for total scattering energy at 2.5 eV. The symbols o and 0 
represent the results obtained by neglecting nuclear spin and 
using Bose-Einstein statistics for the nuclear spin, respectively. 
The discrete data points are connected by a straight line for 
clarity. 

Fig. 2 shows the cross section o v y^ V j for scattering from 
an (v = 0,7 = 0) initial ro-vibrational state to various final 
ro-vibrational states. The curves of o 0 j^oo for low-lying 
final ro-vibrational states, such as (v' = 0,7' = 0) and (v / = 
0,7' = 2) states, have weak oscillating behavior arising 
from the formation of short-lived metastable N 3 states 
[10,1 1]. Note that the cross section for j = 0 f — 48 is 
larger than for j = 0 j' =2 at most energies. The 
variation of the large cross sections for producing high j' 
with kinetic energy verifies that the N + N 2 exchange 
reaction is quite efficient for producing high rotational 
states of N 2 molecules. 

The kinetic energies and average energies of the rotational 
state distributions of our studies are comparable to the 
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translational and rotational temperatures of the 
experiment of [1]. For example, in our study the average 
N 2 rotational energy for a total scattering energy at 3.0 eV 
corresponds to rotational temperature ( T r ) at 5,200 K; the 
kinetic energy in Fig. 2 corresponds to translational 
temperatures in the range 33,000 K to 50,000 K. 
Considering that the vibrational and rotational motions 
are not coupled for low vibrational excitations, we 
conclude that T r for the lower vibrational states is close to 
the value 4,500 K found from the rotational distribution in 
[1]. Furthermore, the translational temperature of the 
experiment, estimated [1] to be above 40,000 K, falls 
within our estimate of 33,000 K to 50,000 K. 




FIG. 3. Cross sections for exchange scattering with vibrational 
excitation from initial ground rotational and vibrational state of 
N 2 , N + N 2 (v = 0, j = 0) — > N 2 (v') + N as a function of kinetic 
energy; v' = 0-3. 

3. C 2 H + H 2 REACTION 

Hydrogen is a major component in the atmosphere of all 
giant planets and methane is a minor component. The 
Titan atmosphere also contains 2-8 percent methane. The 
ethynyl radical C 2 H is an intermediate product of methane 
chemistry upon high temperature shock. The H 2 + C 2 H 
reaction is an important step leading to acetylene, C 2 H 2 . 


FIG. 2. State-to-state cross sections for exchange scattering with 
rotational state excitation from initial ground rotational and 
vibrational state of N 2 , N + N 2 (v =0, j = 0) N 2 (v’ = 0,/) + 

N, vs./ for total scattering energy at 2.5 eV. 

The above results for rotational excitation indicate that the 
exchange reaction of N + N 2 collisions provides a 
mechanism to pump N 2 molecules at low rotational states 
into N 2 molecules at higher rotational states. Because this 
exchange reaction has a propensity for producing highly 
excited rotational states, it provides a route that may lead 
to a nonequilibrium rotational temperature behind a 
strong shock. 

2.2 Vibrational Excitation 

The cross section for v' excitation from a (v = 0,y = 0) 
initial ro-vibrational state is obtained by summing over all 
the final rotational states; i.e., 

= j®v'j'<r- 00 * ( 2 ) 

The values of aV«-o o are shown in Fig. 3 for v' = 0-3. 
This figure shows that significant vibrational excitation of 
N 2 is produced by exchange reactions at higher kinetic 
energies. 


H 2 + C 2 H C 2 H 2 + H. (3) 

Because this reaction also plays an important role in 
combustion chemistry, extensive experimental studies 
have been reported in the literature [12-21], making it a 
good candidate for benchmarking the capability of first 
principles calculations of radical reactions of interest to 
giant planet missions. For this study, the first ever seven- 
dimensional (7D) quantum dynamics calculation is 
carried out. This is done by fixing two of the nine 
dimensional space based on quantum chemistry transition 
state calculation, thus reducing the nine-degree-of- 
freedom calculation to a seven-degree-of-freedom 
calculation. Then, using the time-dependent wave packet 
method [22,23], we calculate the initial state selected 
reaction probabilities for angular momentum J= 0; these 
results are used with /-shifting and energy averaging to 
determine the rate constant. Fig. 4 compares the thermal 
rate constants from the present calculation with those 
from various experimental measurements. It is seen that 
in most cases the calculated rate constants not only agrees 
with the measured data to within experimental error, but 
also reproduce the temperature dependence over a wide 
range of temperatures. There are two exceptions. The 
single measurement of Renlund et al. [13] disagrees with 
both the calculated rate constant and all other 
measurements in the temperature range by almost an 
order of magnitude. Also, the high-temperature data of 
Kruse and Roth [19] not only are higher than the 
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calculated curve but also are higher than the general 
temperature dependence extrapolated from the lower 
temperature data. The present result clearly demonstrates 
our 7D reduced-dimensional calculations can yield 
accurate rate constants for radical reactions. 
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In general, calculations were carried out in two steps. 
First, we used a simulated annealing procedure based on a 
MD method. For complete equilibration, the temperature 
of the system was first increased and then linearly cooled 
down to a low temperature limit. Typical runs in these 
cases were about a picosecond long. The final 
configuration of the low- temperature result was then used 
as an input to a minimization routine based on a quasi- 
Newton algorithm to find the local minimum. This 
method is very efficient in finding low-lying minima for 
systems containing large numbers of moving particles. 

Fig. 5 illustrates the mechanisms of surface defect etching 
of Si(001)-2xl surface using MD simulations. Our results 
show that MD calculations provide an atomic level 
understanding of the etching process, which is the 
removal of halogenated Si particles from the surface 
region. Furthermore, the outcome of these studies 
strongly indicates that MD simulations are able to predict 
the thermodynamics and energetics of surface etching. 


FIG. 4. Calculated rate constant of H 2 + C 2 H C 2 H 2 + H vs. 
inverse temperature. Also presented are various experimental 
results from 0 [12], ♦ [13], x [14], * [15], > [16], V [17], < [18], 
A [19], □ [20], and o[21]. 


4. ETCHING OF A Si SURFACE 

The etching of a Si surface by halogens is used to 
illustrate our approach to the study of surface catalytic 
reactions. In these studies molecular dynamics (MD) 
simulations were used to investigate etching of Si surfaces 
by halogen species [24,25]. For complete descriptions of 
surface reactions taking place during the etching (and 
epitaxial growth) process, it is necessary to understand the 
bonding structures and stabilities of chemisorbed halogen 
species on the Si surface. MD studies today can reproduce 
short-time behavior of an etching process quite 
successfully. Generally speaking, simulation calculations 
are consistent with experiments. 

MD calculations provide detailed understanding of the 
local physics and chemistry taking place during the 
adsorption and subsequent desorption of halogen species 
from the exposed Si surfaces. First and second 
chlorination (and fluorination) steps of surface Si atoms 
were analyzed separately, and the role played by surface 
vacancies was investigated. Steric effects coming from 
repulsive forces among the halogen atoms were found to 
be quite significant. Results indicate that the second step 
that leads to the formation of doubly chlorinated species 
is very important in the overall etching process. In MD 
simulations, energies and forces are calculated employing 
existing empirical potential functions for systems 
containing Si and halogen species. These potentials are 
based on two- and three-body interactions. 


(a) 



FIG. 5. Etching of SiF 3 from A type step edge, (a) The side and 
front views of the sample at the beginning of a room 
temperature MD trajectory calculation, (b) Same as in (a) at the 
end of the MD trajectory. The arrows indicate the dimerization 
of the lower step edge Si atoms as the SiF 3 on the top is etched 
away. 

5. DYNAMIC DATABASE FOR DIATOMICS 

A proper account of re-entry heating requires an accurate 
description of radiation effects. The radiation field is due 
to emission from atoms and molecules that are either 
heated by the shock layer or formed by reactions. Both 
neutrals and ions will participate. The transitions of 
importance are electronic, and modeling requires a 
database of opacity data for each species expected to be 
present in the flow field. 
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Our present focus is on the opacity data of diatomics. The 
diatomics involved depend on the entry atmosphere, and 
the number of species can be quite large. For example, 
entry into the earth's atmosphere produces excited states 
of CN, NO, N 2 , as well as their ions. 

The conventional way to obtain the opacity data is from 
experimental analysis. This, however, is a slow process 
requiring many man-years for a single diatomic. The 
measurements produce line positions and line strengths. 
While the line positions can be obtained with a very high 
degree of accuracy, line strengths are much more difficult 
to measure accurately. In fact, many spectroscopic studies 
do not even attempt to measure line strengths. 

Another source of opacity data is from calculations 
carried out from first principles. In contrast to 
experimental measurements, theory can yield accurate 
line strengths without too much difficulty, but accurate 
line positions are much more difficult. However, progress 
has been made on that front; see, for example, Polyansky 
et al. [26]. In the past, a calculation for a single diatomic 
also required a number of man-years of work. Now we 
can exploit the recent increases in computer power to 
introduce a new paradigm whereby automatic 
computation of accurate opacities for diatomics becomes 
possible. We call our software package for carrying this 
out DDD. 

DDD will be accessible as a world-wide web interface, 
and will work in one of two ways. The first mode, the 
retrieval mode, occurs when a user requests data for a 
diatomic for which calculations had been previously 
carried out. In this case, DDD will return the data 
essentially instantly and in one of several forms indicated 
by the user. The second mode, the data creation mode, 
occurs when a diatomic with no stored data is requested. 
In this case DDD will start the necessary calculations to 
determine the opacity required by the user. Depending on 
the diatomic and the load on the parallel computer used, 
this calculation may require a few minutes to several 
days. If the calculations are not completed in a short 
period of time, DDD will inform the user that calculations 
are being carried out, and will e-mail the user when the 
calculations are completed. 

We have been testing our preliminary version of DDD, 
and the results are quite promising. In Table I we present 
some preliminary results for CN, T e , the electronic 
excitation energy, co e , the vibrational frequency in the 
harmonic approximation, and B e , the rotational constant 
for a rigid rotor. These are low-order spectroscopic 
constants describing the line positions, but our 
calculations also yield line strengths. Numbers in 
parenthesis are experimental values from Huber and 
Herzberg [27]. The agreement is quite satisfactory for this 
level of calculation. Of interest is the large number of 
states that are predicted for which experiments have not 
been carried out. Note that these calculations were carried 


out without any user intervention after the initial 
specification of the atoms C and N. 

The algorithm for DDD is quite complicated. Given a 
database for atoms consisting of orbital occupations of the 
ground and valence excited states, DDD constructs a 
preliminary list of binding configurations of the diatomic. 
These orbital occupations are then used in trial 
calculations. Based on these results, more sophisticated 
calculations are carried out on the strongly bound levels 
— potential energy curves, transition moment curves, and 
coupling matrix elements are computed. For these 
calculations, the highly reliable multi-reference 
configuration interaction method is used. These are then 
used in a coupled rotation- vibrational-electronic code [28] 
to generate a line list. 


TABLE I: Summary of Electronic States of CN. Numbers 
in parenthesis are from [27]. 



U cm' 1 ) 

0) e 

B e 

4 2 + 

67917 

1073 

1.316 

2 n r 

67620 

1351 

1.174 

4 n r 

66585 

2393 

1.714 

2 A r 

60704 (60096) 

1549(1239) 

1.375(1.383) 

2 r 

59638 

1242 

1.347 

2 s + 

59439(59151) 

1721 (1681) 

1.456(1.487) 

2 O r 

59291 

931 

1.095 

2 n r 

58747 

2187 

1.429 

4 r 

52655 

1260 

1.348 

2 n r 

51761 (54486) 

1133 (1004) 

1.045(1.162) 

% 

46864 

1309 

1.367 

4 n r 

44243 

923 

1.104 

4 2 + 

36466 

1326 

1.367 

2 2 + 

27128 (25752) 

2059 (2163) 

1.868(1.973) 

2 n ; 

8401 (9245) 

1766(1812) 

1.642(1.715) 

2 s + 

0(0) 

2028 (2068) 

1.818(1.900) 


6. CONCLUDING REMARKS 


This paper illustrates the capability of present day first 
principles quantum chemistry calculations. By making 
use of the recent advances in hardware and software, first 
principles calculations can now provide reliable chemistry 
data needed for aerothermal analysis of aerocapture 
missions and other entry problems. The N + N 2 
calculation supports the experimental observation by 
Fujita et al. [1] of strong nonequilibrium between the 
kinetic and rotational temperatures of a high-energy 
shock. Furthermore, our theoretical calculations provide 
insight into the origin of this nonequilibrium and indicate 
that future aerothermal modeling of high-velocity entry 
should incorporate these nonequilibrium conditions. The 
scatter in the measured C 2 H + H 2 reaction data at high 
temperatures illustrates the risk of using limited 
experimental data. First principles calculations can be 
validated by comparisons with available low- temperature 
experimental data and, moreover, provide otherwise 
unavailable high-temperature data. Our capability is not 
limited to gas-phase studies, as illustrated by the silicon 
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etching study. These techniques can be used to study 
surface catalytic reactions. The web-based DDD is a 
pioneering tool that provides a powerful means for the 
aerothermodynamicists to do their modeling of radiation 
heat load. We also point out that sensitivity analysis 
should be used as a guide for selecting the critical subset 
of required data. This approach will significantly shorten 
the turn around time, minimize the computational effort, 
and maximize the benefit of the calculations. 
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